Identification of the dominant precession damping mechanism in Fe, Co, and Ni by 

first-principles calculations 
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The Landau-Lifshitz equation reliably describes magnetization dynamics using a phenomenolog- 
ical treatment of damping. This paper presents first-principles calculations of the damping param- 
eters for Fe, Co, and Ni that quantitatively agree with existing ferromagnetic resonance measure- 
ments. This agreement establishes the dominant damping mechanism for these systems and takes 
a significant step toward predicting and tailoring the damping constants of new materials. 



Magnetic damping determines the performance of 
magnetic devices including hard drives, magnetic ran- 
dom access memories, magnetic logic devices, and mag- 
netic field sensors. The behavior of these devices can be 
modeled using the Landau-Lifshitz (LL) equation [l| 



m = -|7|m X Hcff 



X (m X Heff) , (1) 



or the essentially equivalent Gilbert (LLC) form 0, 
The first term describes precession of the magnetization 
m about the effective field Hoff where 7 = is the 

gyromagnetic ratio. The second term is a phenomeno- 
logical treatment of damping with the adjustable rate A. 
The LL(G) equation adequately describes dynamics mea- 
sured by techniques as varied as ferromagnetic resonance 
(FMR) [4] , magneto-optical Kerr effect 0] , x-ray absorp- 
tion spectroscopy 6], and spin-current driven rotation 
with the addition of a spin-torque term 0, Q . 

Access to a range of damping rates in metallic mate- 
rials is desirable when constructing devices for different 
applications. Ideally, one would like the ability to de- 
sign materials with any desired damping rate. Empiri- 
cally, domng NiFe alloys with transition metals ;9,] or rare 
earths [10| has produced compounds with damping rates 
in the range of a = 0.01 to 0.8. A recent investigation 
of adding vanadium to iron resulted in an alloy with a 
damping rate slightly lower than that for pure iron llj, 
the system with the lowest previously known value. How- 
ever, the damping rate of a new material cannot be pre- 
dicted because there has not yet been a first-principles 
calculation of damping that quantitatively agrees with 
experiment. The challenging pursuit of new materials 
with specific or lowered damping rates is further com- 
plicated by the expectation that, as device size contin- 
ues to be scaled_down, material parameters, such as A, 
should change [l^. A detailed understanding of the im- 
portant damping mechanisms in metallic ferromagnets 
and the ability to predictively calculate damping rates 
would greatly facilitate the design of new materials ap- 
propriate for a variety of applications. 

The temperature dependence of damping in the tran- 
sition metals has been carefully characterized through 
measurement of small angle dynamics by FMR [l3|. 



While one might naively expect damping to increase 
monotonically with temperature, as it does for Fe, both 
Co and Ni also exhibit a dramatic rise in damping at low 
temperature as the temperature decreases. These ob- 
servations indicate that two primary mechanisms are in- 
volved. Subsequent experiments [ij, [3] partition these 
non-monotonic damping curves into a conductivity-like 
term that decreases with temperature and a resistivity- 
like term that increases with temperature. The two terms 
were found to give nearly equal weight to the damping 
curve of Ni and have temperature dependencies similar 
to those of the conductivity and resistivity, suggesting 
two distinct roles for electron-lattice scattering. 

The torque-correlation model of Kambersky [31 ap- 
pears to qualitatively match the data. However, like 
most of the various models presented by Kambersky 
[3 [13, E, Hi and others [20|, it has not been quan- 
titatively evaluated in a rigorous fashion. This has left 
the community to speculate, based on rough estimates or 
less, as to which damping mechanisms are important. We 
resolve this matter in the present work by reporting first- 
principles calculations of the Landau-Lifshitz damping 
constant according to Kambersky's torque-correlation 
expression. Quantitative comparison of thejDresent cal- 
culations to the measured FMR values 13] positively 
identifies this damping pathway as the dominant effect 
in the transition metal systems. In addition to present- 
ing these primary conclusions, we also describe the re- 
lationship between the torque-correlation model and the 
more widely understood breathing Fermi surface model 
[2l[, showing that the results of both models agree 
quantitatively in the low scattering rate limit. 

The breathing Fermi surface model of Kambersky pre- 
dicts 
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This model offers a qualitative explanation for the low 
temperature conductivity-like contribution to the mea- 
sured damping. The model describes damping of uniform 
precession as due to variations den,k/d9 in the energies 
en,fc of the single-particle states with respect to the spin 
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direction 9. The states are labeled with a wavevector 
k and band index n. As the magnetization precesses, 
the spin-orbit interaction changes the energy of electronic 
states pushing some occupied states above the Fermi level 
and some unoccupied states below the Fermi level. Thus, 
electron-hole pairs are generated near the Fermi level 
even in the absence of changes in the electronic popula- 
tions. The T] function in Eq. ([2]) is the negative derivative 
of the Fermi function and picks out only states near the 
Fermi level to contribute to the damping, g is the Lande 
g-factor and /zb is the Bohr magneton. The electron-hole 
pairs created by the precession exist for some lifetime r 
before relaxing through lattice scattering. The amount of 
energy and angular momentum dissipated to the lattice 
depends on how far from equilibrium the system gets, 
thus damping by this mechanism increases linearly with 
the electron lifetime as seen in EqlS] Since the electron 
lifetime is expected to decrease as the temperature in- 
creases, this model predicts that damping diminishes as 
the temperature is raised. 

Because the predicted damping rate is linear in the 
scattering time the damping rate cannot be calculated 
more accurately than the scattering time is known. For 
this reason it is not possible to make quantitative com- 
parisons between calculations of the breathing Fermi sur- 
face and measurements. Further, while the breathing 
Fermi surface model can explain the dramatic temper- 
ature dependence observed in the conductivity-like por- 
tion of the data it fails to capture the physics driving the 
resistivity-like term. This is a significant limitation from 
a practical perspective because the resistivity-like term 
dominates damping at room temperature and above and 
is the only contribution observed in iron and NiFe 
alloys (l^l. For these reasons it is necessary to turn to 
more complete models of damping. 

Kambersky's torque-correlation model predicts 

^ = ^E / |r_(fc)|V_(fc) (3) 

and we will show that it both incorporates the physics of 
the breathing Fermi surface model and also accounts for 
the resistivity-like terms. The matrix elements r~^(fc) = 
(n, k\[a~ , Hso]\m, k) measure transitions between states 
in bands n and m induced by the spin-orbit torque. 
These transitions conserve wavevector k because they de- 
scribe the annihilation of a uniform precession magnon, 
which carries no linear momentum. The nature of these 
scattering events, which are weighted by the spectral 
overlap = (I/tt) / dwi 7?(wi)A„fc(a;i)Amfe(a;i), 

will be discussed in more detail below. The electron spec- 
tral functions Ank are Lorentzians centered around the 
band energies e„fc and broadened by interactions with the 
lattice. The width of the spectral function Tl/t provides a 
phenomenological account for the role of electron-lattice 
scattering in the damping process. The 77 function is the 
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FIG. 1: Calculated Landau-Lifshitz damping constant for Fe, 
Co, and Ni. Thick solid curves give the total damping param- 
eter while dotted curves give the intraband and dashed lines 
the interband contributions. The top axis is the full-width- 
half-maximum of the electron spectral functions. 



same as in Eq. ^ and enforces the requirement of spec- 
tral overlap at the Fermi level. 

Equation ^ captures two different types of scatter- 
ing events: scattering within a single band, m = n, for 
which the initial and final states are the same, and scat- 
tering between two different bands, m^n. As explained 
in |16j the overlap of the spectral functions is propor- 
tional (inverse) to the electron scattering time for intra- 
band (interband) scattering. From this observation the 
qualitative conclusion is made that the intraband contri- 
butions match the conductivity-like terms while the inter- 
band contributions give the resistivity-like terms. While 
this seems promising, evaluation of Eq. Q is more com- 
putationally intensive than that of the breathing Fermi 
surface model and until now only a few estimates for Ni 
and Fe have been made Il9|. 



3 



TABLE I: Calculated and measured [TH damping parameters. Values for A, the Landau-Lifshitz form, are reported in 10^ , 
values of a, the Gilbert form, are dimensionless. The last two columns list calculated damping due to the intraband contribution 
from Eq. Q and from the breathing Fermi surface model respectively. Values for A/r are given in 10^^ s-^ Published 
numbers from [l3l ] and have been multiplied by 4n to convert from the cgs unit system to SI. 
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bcc Fe (001) 


0.0013 


0.54 
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0.61 
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0.968 


bcc Fe (111) 


0.0013 


0.54 






1.35 


1.29 


hep Co (0001) 


0.0011 


0.37 


0.9 


0.41 


0.786 


0.704 


fee Ni (111) 


0.017 


2.1 


2.9 


0.72 


6.67 


6.66 


fee Ni (001) 


0.018 


2.2 






8.61 


8.42 



We have performed first-principles calculations of the 
torque-correlation model Eq. ([3]) with realistic band 
structures for Fe, Co, and Ni. Prior to evaluating Eq. ([3]) 
the eigenstates and energies of each metal were found us- 
ing the linear augmented plane wave method [ 2 3| in the 
local spin density approximation (LSDA) [24l. |25|. [2^. 
Details of the calculations for these materials are de- 
scribed in [23] ■ The exchange field was fixed in the cho- 
sen equilibrium magnetization direction. Calculations of 
Eq. ([3]) presented in this paper are converged to within 
a standard deviation of 3 %, which required sampling 
(160)3 fc-points for Fe, (120)^ for Ni, and (100)^ fc-points 
in the basal plane by 57 along the c-axis for Co. Electron- 
lattice interactions were treated phenomenologically as a 
broadening of the spectral functions. The Fermi distribu- 
tion was smeared with an artificial temperature. Results 
did not vary significantly with reasonable choices of this 
temperature since the broadening of the Fermi distribu- 
tion was considerably less than that of the bands. The 
damping rate was calculated for a range of scattering 
rates (spectral widths) just as damping has been mea- 
sured over a range of temperatures. 

The results of these calculations are presented in Fig.[T] 
and are decomposed into the intraband and interband 
terms. The downward sloping line in Fig. [1] represents 
the intraband contribution to damping. Damping con- 
stants were rece ntly calculated using the breathing Fermi 
surface model [12 [2l[ by evaluating the derivative of 
the electronic energy with respect to the spin direc- 
tion according to Eq. ([2]). The results of the breathing 
Fermi surface prediction are indistinguishable from the 
intraband terms of the present calculation even though 
the computational approaches differed significantly; the 
agreement is quantified in Table [H 

The breathing Fermi surface model could not be quan- 
titatively compared to the experimental results because 
the temperature dependence of the scattering rate has 
not been determined sufficiently accurately. While the 
present calculations also require knowledge of the scat- 
tering rate to determine the damping rate the non- 



monotonic dependence of damping on the scattering rate 
produces a unique minimum damping rate. In the same 
manner that the calculated curves of Fig. [1] have a mini- 
mum with respect to scattering rate, the measured damp- 
ing curves exhibit minima with respect to temperature. 
Whatever the relation between temperature and scatter- 
ing rate, the calculated minima may be compared di- 
rectly and quantitatively to the measured minima. Ta- 
ble [T] makes this comparison. The agreement between 
measured and calculated values shows that the torque- 
correlation model accounts for the dominant contribution 
to damping in these systems. 

Our calculated values are smaller than the measured 
values. Using measured g values instead of setting g — 
2 would increase our results by a factor of (g/2)^, or 
about 10 % for Fe and 20 % for Co and Ni. Other pos- 
sible reasons for the difference include a simplified treat- 
ment of electron-lattice scattering in which the scattering 
rates for all states were assumed equal, the mean-field 
approximation for the exchange interaction, errors asso- 
ciated with the local spin density approximation (LSDA) , 
and numerical convergence (discussed below). Other 
damping mechanisms may also make small contributions 



Since the manipulations involved with the equation of 
motion techniques employed in deriving Eq. ^ obscure 
the underlying physics we now discuss the two scatter- 
ing processes and connect the intraband terms to the 
breathing Fermi surface model. The intraband terms in 
Eq. ([3]) describe scattering from one state to itself by 
the torque operator, which is similar to a spin-flip oper- 
ator. A spin-flip operation between some state and itself 
is only non-zero because the spin-orbit interaction mixes 
small amounts of the opposite spin direction into each 
state. Since the initial and final states are the same, the 
operation is naturally spin conserving. The matrix ele- 
ments do not describe a real transition, but rather pro- 
vide a measure of the energy of the electron-hole pairs 
that are generated as the spin direction changes. The 
electron-hole pairs are subsequently annihilated by a real 
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electron-lattice scattering event. 

To connect the derivatives de/d9 in Eq. ([2]) and the 
torque matrix elements in Eq. ([3]) we imagine first point- 
ing the magnetization in some direction z. The only 
energy that changes with the magnetization direction is 
the spin-orbit energy Hso- As the spin of a single parti- 
cle state I) rotates along 6 about x its spin-orbit energy 
is given by e{e) = (|e*'"="^ i?so e"*'""=''|). The derivative 
with respect to 6 is de{9)/d9 i{\e"'-''[a^ , i?so]e~*'"^^|). 
Evaluating this derivative at the pole {9 = 0) gives 
de/d9 = i{\[ax , Hso]\)- Similarly, rotating the spin along 
9 about y leads to de/dO — i{\[(Jy , -ffsolD- The torque 
matrix elements in Eq. ^ are F" = (|[a-~ , i?so]|) = 
{\[<Jx , Hso\\) —i{\[<7y , -f^solD- Using the relations between 
the commutators and derivatives just found the torque is 
= —i[de/d9)x — {de/d9)y where the subscripts in- 
dicate the rotation axis. Squaring the torque matrix 
elements gives \T-\^ = {de/d9)l + {de/dfffy. For high 
symmetry directions {de/d9)x — {de/d9)y and we de- 
duce |F~p — 2{de/d9Y demonstrating that the intra- 
band terms of the torque-correlation model describe the 
same physics as the breathing Fermi surface. 

The monotonically increasing curves in Fig. [1] indi- 
cate the interband contribution to damping. Uniform 
mode magnons, which have negligible energy, may in- 
duce quasi-elastic transitions between states with differ- 
ent energies. This occurs when lattice scattering broad- 
ens bands sufficiently so that they overlap at the Fermi 
level. These wavevector conserving transitions, which are 
driven by the precessing exchange field, occur primarily 
between states with significantly different spin character. 
The process may roughly be thought of as the decay of 
a uniform precession magnon into a single electron spin- 
fiip excitation. These events occur more frequently as 
the band overlaps increase. For this reason the interband 
terms, which qualitatively match the resistivity-like con- 
tributions in the experimental data, dominate damping 
at room temperature and above. 

We have calculated the Landau-Lifshitz damping pa- 
rameter for the itinerant ferromagnets Fe, Co, and Ni as a 
function of the electron-lattice scattering rate. The intra- 
band and interband components match qualitatively to 
conductivity- and resistivity-like terms observed in FMR 
measurements. A quantitative comparison was made be- 
tween the minimal damping rates calculated as a function 
of scattering rate and measured with respect to temper- 
ature. This comparison demonstrates that our calcula- 
tions account for the dominant contribution to damping 
in these systems and identify the primary damping mech- 
anism. At room temperature and above damping occurs 
overwhelmingly through the interband transitions. The 
contribution of these terms depends in part on the band 
gap spectrum around the Fermi level, which could be 
adjusted through doping. 
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